Interplay between geometry and temperature for inclined Casimir plates 
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We provide further evidence for the nontrivial interplay between geometry and temperature in the 
Casimir effect. We investigate the temperature dependence of the Casimir force between an inclined 
semi-infinite plate above an infinite plate in D dimensions using the worldline formalism. Whereas 
the high-temperature behavior is always found to be linear in T in accordance with dimensional- 
reduction arguments, different power-law behaviors at small temperatures emerge. Unlike the case 
of infinite parallel plates, which shows the well-known behavior of the force, we find a T^^^ 
behavior for inclined plates, and a ~ ^gji^vior for the edge effect in the limit where the plates 

become parallel. The strongest temperature dependence ~ y^-^ Q^^^^g f^j. ^j^g Casimir torque of 
inclined plates. Numerical as well as analytical worldline results are presented. 



I. INTRODUCTION 

The Casimir effect [l| is not only a field witnessing 
rapid experimental as well as theoretical progress, it also 
continues to offer surprising new features. The Casimir 
effect derives its fascination from the fact that it origi- 
nates from quantum fluctuations of the radiation field or 
of the charge distribution on the mesoscopic or macro- 
scopic test bodies. Moreover, it inspires many branches 
of physics, ranging from mathematical to applied physics, 
see [2,S| for reviews and [1] for experimental verifications. 

A distinctive feature of Casimir forces between test 
bodies is the dependence on the geometry, i.e., the shape 
and orientation of these bodies. For a comparison be- 
tween theory and a real Casimir experiment, a number 
of properties such as finite conductivity, surface rough- 
ness and finite temperature have to be taken into account 
in addition. Generically, these latter corrections do not 
factorize but take infiuence on one another. For instance, 
the interplay between dielectric material properties and 
finite temperature [3| is still a subject of intense theoret- 
ical investigations and has created a long-standing con- 
troversy 0,0, H, Q • Also the role of electrostatic patch 
potentials has been suggested as a potentially problem- 
atic iss ue flol, [l]] |. which has become a matter of severe 
debate 

The present article is not meant to resolve these con- 
troversies. On the contrary, our work intends to draw 
attention to another highly nontrivial interplay which 
on the one hand needs to be accounted for when com- 
paring theory and a real experiment and on the other 
hand is another characteristic feature of the Casimir ef- 
fect: the interplay between geometry and temperature. 
As first conjectured by Jaffe and Scardicchio the 
temperature dependence of the Casimir effect can be 
qualitatively different for different geometries, as both 
the pure Casimir effect as well as its thermal corrections 
arise from the underlying spectral properties of the fiuc- 
tuations. First analytical as well as numerical evidence 
of this "geothermal" interplay in a perpendicular-plates 
configuration has been found in using the worldline 



formalism. 

The physical reason for this interplay can be under- 
stood in simple terms: for the classical parallel-plate case, 
the nontrivial part of the fiuctuation spectrum is given 
by the modes orthogonal to the plates. This relevant part 
of the spectrum has a gap of WcLve number fcgap 

= ir/a, 

where a is the plate separation. For small temperatures 
^gap) the higher- lying relevant modes can hardly 
be excited, such that their thermal contribution to the 
Casimir force remains suppressed: the resulting force law 
for the parallel-plates case scales like (aT)^. This argu- 
ment for a suppression of thermal contributions applies 
to all geometries with a gap in the relevant part of the 
spectrum (e.g. concentric cylinders or spheres, Casimir 
pistons, etc.). These geometries are called closed.^ 

This reason for a suppression of thermal contributions 
is clearly absent for open geometries with a relevant gap- 
less part of the spectrum. For these geometries, rele- 
vant modes of the spectrum can always be excited at any 
small temperature value. Therefore, a stronger thermal 
contribution ~ (aT)°' with < a < 4 can be expected. 
As experimentally important configurations such as the 
sphere-plate or the cylinder-plate geometry belong to this 
class of open geometries, a potentially significant geother- 
mal interplay may exist in the relevant parameter range 
aT- 0.01... 0.1. 

In the present work, we provide further evidence for 
the geometry-temperature interplay in the Casimir effect. 
For simplicity, we study the Casimir effect induced by a 
fiuctuating real scalar field obeying Dirichlet boundary 
conditions ('Dirichlet scalar'). As an illustrative exam- 
ple, we concentrate on an inclined-plates configuration; 
here, a semi-infinite plate is located above an infinite one, 
with an angle of inclination of < < 7r/2, see Fig. [TJ 



^ Of course, parallel plates as well as concentric cylinders are not 
closed in the sense of compactness. Also, they have a gapless 
part of the spectrum along the symmetry axes. However, this 
part of the spectrum does not give rise to the Casimir force and 
hence is not a relevant part. 
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gral approach |39l.l40l.l4ll|. It will certainly be worthwhile 
to generalize these methods to finite temperature for a 
study of the geometry-temperature interplay. 

In the remainder of this introduction, we summarize 
our most important results specializing to 3 + 1 dimen- 
sional spacetime. In Sect-Ull we briefiy review the world- 
line approach to the Casimir effect. Sect. [Till is devoted 
to a study of the zero-temperature Casimir effect for the 
geometries under consideration. The finite-temperature 
case is described in Sect. HVl Our conclusions are sum- 
marized in Sect. IVl 



A. Summary of results in D = 4 



FIG. 1: Sketch of the inclined-plates configuration. The infi- 
nite plate (dashed line) is rotated in the x, z plane by an angle 
f). As special cases, — ^ corresponds to the configuration 
of one semi-infinite plate parallel to an infinite plate (Isi con- 
figuration), whereas = 7r/2 yields the perpendicular-plates 
configuration. 



This configuration generalizes geometries which have first 
been proposed and studied in the context of Casimir edge 
effects Our results do not only generalize the find- 
ings of jl5| which hold for — it/2. Most importantly, we 
identify regimes with fractional temperature dependences 
for certain geometries. Moreover, we work m D = d + 1 
dimensional spacetime, yielding many analytical as well 
as numerical results for the Casimir force and energy as 
well as for the torque. 

A reliable study of geothermal Casimir phenomena re- 
quires a method that is capable of dealing with very 
general Casimir geometries. For this, we use the world- 
line approach to the Casimir effect [17||, which is based 
on a mapping of field-theoretic fiuctuation averages onto 
quantum- mechanical path integrals (isl. fiol. [20| . For ar- 
bitrary backgrounds, this worldline integral represent- 
ing the spacetime trajectories of the quantum fiuctua- 
tions can straightforwardly be computed by Monte Carlo 
methods [2l|. As the computational algorithm is gen- 
erally independent of the background, i.e., the Casimir 
geometry in our case, Casimir problems can straightfor- 
wardly be tackled with this method. High-precision com- 
putations for Dirichlet-scalar fiuctuations have been per- 
formed, e^., for the sphere-plate and cylinder-plate case 
[13, [13, 13 ■ III the present work, we demonstrate that 
the worldHne approach can also be used to obtain novel 
analytical results (see also for an analytical worldline 
approximation technique) . 

In order to overcome standard approximative tools 
based, e.g., on the proximity-force theorem [22], a va- 
riety of new field-theoretical methods for Casimir phe- 
nomena have been developed in recent yearSj rang- 
ing from improved approximation methods [13, El, 0| 
to exact methods mainly based on scattering theory 
S [U, [13, [13, [13, [H, m, [13, [Hi or a functional inte- 



Let us already summarize our most important results, 
specializing to 13 = 4 spacetime dimensions and con- 
centrating on the Casimir interaction energy; the corre- 
sponding force can straightforwardly be derived by dif- 
ferentiation. At zero temperature, the classical Casimir 
energy of two parallel Dirichlet plates at a distance a 
reads 



A 



cuhc 



1440 



0.00685, 



(1) 



where A is the area of the plates. From now on, we use 
natural units, setting he = I. The Casimir energy of 
inclined plates (i.p.) can be parameterized as 





L,, 



sm{ip) a? ' 



(2) 



where the coefficient c<^ is shown in Fig. |4] as a function 
of (/3. The extent of the inclined plate in y direction along 
the edge is Ly. At (/? = 0, the energy per edge length ^ 
diverges and has to be replaced by 



E 



E. 



lsi,cdgc 



A 



EyCQdg 



(3) 



where Ec''^ is the Casimir Energy |T| with A^^^ being 
the semi-infinite plate's area and E^^^'"^^'^ the so-called 
edge energy. The numerical value of Codgc is about 0.0026 
in agreement with (l6| . 

The Casimir torque is obtained from Eq. ^ by 
£)i.p.,v = dS^P-'^/d^. For £> = 4, the torque D'^P-'f as a 
function of (p is shown in Fig. [H At </? = 0, the Casimir 
torque per unit length diverges as well but can be con- 
verted into finite torque per unit area. Remarkably, for 
(f = the standard torque obtained from the Casimir 
energy of parallel plates ^ dI = AL^-k^ /960a-^ « 
0.0103^Lz/a* is reduced by a repulsive contribution 
« — 0.003660iy/a^ arising from the edge effect. We en- 
counter a similar subleading repulsive torque effect at 
finite temperature. 

Thermal fiuctuations modify the Casimir energy, yield- 
ing the free energy 



E,iT)^Ec{0) + AE^{T), 



(4) 
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where AEc{T) is the temperature correction. For (aT) 
0, the correction AEc {aT to the well-known parallel- 
plates energy reads 



AE!{aT^ 0) 



47r 



90 



(5) 



which is « -0.0957T3 + 0.110ar4. Note that only the 
term contributes to the force as the first term vanishes 
upon differentiation. 

For (aT) 0, our result for the thermal correction 
Ai?^ P-''^(T) to the inclined-plates energy reads 



A£;^-P-'^(aT^ 0) 



C(3)ar3 
24sin(iy9) ' 47rsin((^)' 



(6) 



where Ci^^To shown in Fig. [9] as a function of (p. The 
second term which is a purely analytical result is the 
generalization of a result for perpendicular plates, ip = 
7r/2, found in numerically, this term evaluates to 

w 0.0957arVsin((/9). 

Again, Eq. ^ denotes an energy per edge length and 
diverges as </? ^ 0. It has to be replaced by the formula 
for the energy of a semi-infinite plate above a parallel 
one, El'^'iT) = £;i>5i^<=dgo^y) + i;c'''"(T). The thermal 
part of £'c'^''"(T) is as in ([5]), where A is the area of the 
semi-infinite plate. The leading thermal correction to the 
edge effect AEl^'''"^s°{T) reads 



^^Isi.odgoj-ji-j 



-^^|^+0.063al•^4T3■^^ (7) 



For [aT] —* oo, all thermal Casimir energies increase 
linearly in T due to dimensional reduction. For instance. 



the Casimir energy Ec (T) for parallel plates becomes 



eI {aT ^ oo) = - 



AC{i)T 



Sna^ 



(8) 



which is « -0M78AT/a^. Note that Ej^iaT oo) is in- 
dependent of he as the dimensional analysis easily shows. 
The energy at large (aT) can therefore be interpreted as 
a classical effect. 

The same holds for the large (aT) behavior of the 
inclined-plates case as well as for semi-infinite plates. For 
inclined plates, we get 

El:^-'^iaT^oo) = -^^f{p^^, (9) 
(47r)2a sm((^) 

where c^p^T^ is shown in Fig. [9] as a function of ip. 
The edge effect reads at large (aT) 

^l.i,cdgc(^y ^ oo) = (10) 

a 

In main part of this article, these results will be derived 
in detail in D spacetime dimensions. 



II. WORLDLINE APPROACH TO THE 

CASIMIR EFFECT 

Let us briefiy review the worldline approach to the 
Casimir effect for a massless Dirichlet scalar; for details, 
see (itI. [23|. Consider a configuration S consisting of two 
rigid objects with surfaces Si and S2. The worldline rep- 
resentation of the Casimir interaction energy in D = d+ 1 
dimensional spacetime reads 



Ec 



1 



dT 



2(47r)^/2 7o 



dVM(es[x(r)]). (11) 



Here, the generalized step functional obeys Gx;[x] = 1 if 
a worldline x(t) intersects both surfaces S = Ei U S2, 
and is zero otherwise. 

The expectation value in Eq. ifTTj) is taken with respect 
to an ensemble of d-dimensional closed worldlines with 
a common center of mass xqm and a Gaufiian velocity 
distribution, 



f Dx 



(12) 



Here, we have already used the fact that the time compo- 
nent cancels out for static Casimir configurations at zero 
temperature. Eq. ifTTj) has an intuitive interpretation: 
All worldHnes intersecting both surfaces do not satisfy 
Dirichlet boundary conditions on both surfaces. They 
are removed from the ensemble of allowed fiuctuations 
by the 8 functional and thus contribute to the negative 
Casimir interaction energy. In the process of the aux- 
ihary T integration, the propertime parameter T scales 
the extent of a worldline by a factor of -y/T. Large T 
correspond to long-wavelength or IR fiuctuations, small 
T to short- wavelength or UV fiuctuations. 

Introducing finite temperature T = 1//3 by the Mat- 
subara formalism is equivalent to compactifying Eu- 
cHdean time on the interval [0, /?]. Now, the closed world- 
lines live on a cylindrical surface and can carry a winding 
number. The worldlines x*^"^ (r) winding n times around 
the cylinder can be decomposed into a worldline x(t) 
with no winding number and a winding motion at con- 
stant speed. 



(13) 



where the Dth component corresponds to EucHdean 
time. The Casimir energy ifTTj) now becomes 



Er 



1 



2(47r)^/2 
dT 



(14) 



7-1+-D/2 



E 



dVivi (es[x(r)]) 



The finite-temperature worldline formalism for static 
configurations thus boils down to a winding-number pref- 
actor in front of the worldline expectation value together 
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with a sum over winding numbers: 



1 



2E^- 

n=l 



(15) 



The winding-number sum is directly related to the stan- 
dard Matsubara sum by a Poisson resummation, 



1 



E 



(16) 



This is already sufficient to understand the high- 
temperature limit of generic Casimir configurations: at 
high temperatures /3 ^ 0, only the zeroth Matsubara fre- 
quency survives as higher modes receive thermal masses 
of order ~ 27r//3 = 27rT and decouple. All remaining tem- 
perature dependence arises from the dimensional prefac- 
tor 1//3 — T, and the dependence on the Casimir ge- 
ometry only enters the prefactor. The calculation of the 
latter is a dimensionally reduced problem in D — 1 di- 
mensions. This is a general mechanism of dimensional 
reduction in high-temperature field theories. The Hnear 
high-temperature asymptotics is also clear from the fact 
that the Bose-Einstein distribution governing the distri- 
bution of bosonic thermal fiuctuations increases as ~ T 
in the high-temperature limit. 

Finally, it is advantageous for numerical as well as an- 
alytical calculations to rescale the worldHnes such that 
the velocity distribution becomes independent of T, 

7(0:=^x(T<) ^g-i/o^'dt^ (17) 

where the dot always denotes a derivative with respect to 
the argument, e.g., 7 = d'y{t)/dt. In terms of these nor- 
malized worldlines 7 and the center-of-mass coordinate 
xcM, the 9 function reads more expHcitly 

e[x] = e[xcM + \/T7(t)]. (18) 

The involved worldline integrals can be evaluated also 
numerically by Monte Carlo methods in a straightforward 
manner. For this, the path integral over an operator O 
is approximated by a sum over a finite ensemble of til 
loops, 



{Ob]) - -E^t-/^]' 



(19) 



where £ counts the worldHnes in the ensemble. Each 
worldline 'y{t) is furthermore discretized by a finite set 
of N points per loop (ppl). 



lit) 



It 



liU), U 



TV' 



0, 



(20) 



where 70 = Jn are identified as the worldlines are closed. 
Various efficient ab initio algorithms for generating dis- 
cretized worldlines with Gaufiian velocity distribution 
have been developed, see, e.g.. 

With these comparatively simple prerequisites, we can 
now turn to an analysis of various non-trivial Casimir 
configurations for the Dirichlet scalar. 
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FIG. 2: Dth moment of the maximum spatial extent A of a 
worldline as a function of D. This geometric object (which is 
the same in any target dimension of the worldline) is related 
to the Casimir energy of the parallel-plates configuration in 
D spacetime dimensions by Eqs. l[23|l and EH The plot com- 
pares the exact analytical result with the worldline numerical 
computation based on 1000 worldlines with 2x 10® ppl (points 
per loop) each. 



III. CASIMIR EFFECT AT ZERO 
TEMPERATURE 

Let us first study parallel and inclined plates at zero 
temperature. The purpose of this section is on the one 
hand to review and generalize known results and on the 
other hand to exemplify how the Casimir effect can be 
understood in terms of simple geometric properties of the 
worldlines. 



A. Parallel Plates 

We start with Casimir's classic configuration of two in- 
finitely extended parallel plates. Let the lower and upper 
plate lie in the z = —a and z = planes, respectively. 
In d space dimensions the surface area A of the plates is 
then d — 1 dimensional. The 8 functional for this config- 
uration reads 



Oil zcM + Vr^z.e = ( VT ^z^^^j + zcM 



(21) 



-ZcM 



riz 



where is the z coordinate of the £'th worldline (mea- 
sured with respect to the center of mass). The quanti- 
ties 7z„a^,f, 7z„,i„,£ denote the worldline's maximal and 
minimal extent in the z direction, respectively. The 
total maximal extent \i of the ^'th worldline then is 
— Iz^^^.i — 7z,„i„,f- Now, we can do the integral in 
Eq. (flTI) and obtain the Casimir energy density (sup- 
pressing the index £ from now on) 



A 



(A 



£)(£)- l)(47r)^/2 



(22) 
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We observe that the D-dimensional parallel-plate 
Casimir energy is related to the Dth cumulant of the 
extent of the worldHnes This is a first example for a 
relation between Casimir energies and geometric proper- 
ties of the worldlines. Instead of computing these cumu- 
lants directly, let us simpl y c ompare Eq. (|22l) with the 
well-known analytic result |43.[44|. 



El 
A 



TiD/2)C{D) 



yielding 



{X^)=D{D-mD/2)aD). 



(23) 



(24) 



A comparison of the analytical result to a numerical eval- 
uation of the cumulants is displayed in Fig. [21 Also, the 
Casimir force density can straightforwardly be obtained 
as the derivative of Eq. I|22p with respect to a. Inciden- 
tally, the connection between Casimir energies and world- 
line properties also induces a relation between Casimir 
energies and questions in polymer physics, as first ob- 
served in (241. 



|7x^inb)l 



T^max (^) 




FIG. 3: All relevant information for the evaluation of the 
Casimir energy of the inclined plates, l[28ll . is encoded in 
the function 7m (a;), which has to be integrated from 7a;„i„ 

to 7a:„,ax ■ 



B. Inclined plates 

The incHned-plates (i.p.) configuration consists of a 
perfectly thin semi-infinite plate above an infinite plate 
at an angle (p, see Fig. [TJ The semi-infinite plate has an 
edge with a (d — 2 dimensional) length Ly. The infinite 
plate has a (d — 1) dimensional area A."^ Let a be the 
minimal distance between the plates. In the following, 
we will omit the center-of-mass subscript CM. The 
functional for this configuration reads 

©i.p..^; =^ i-xcos{ip) - zsm{(p) - VT 7x„i„,f(v) - a) 




where the first function ensures the intersection of the 
worldline with the infinite plate. The remaining three 
ones account for the intersection with the semi-infinite 
plate. In Eq. ((25|) . we have used 

Ix^i^Ay^) = min i^xAi) cos{ip) + ^z,i{t) sin(<^)) , (26) 

where t parameterizes the worldline; i.e., in the dis- 
cretized version, we have t = 1 . . .N with N being the 
number of points per worldline loop (ppl). Trivially, 

7x„i„,f(0) = 7x„i„,£ and 72-„i„/(7i'/2) = j^^.^^i holds. 



^ Of course, the labels "semi-infinite" and "infinite" imply that both 
Ly and A are considered in the limit Ly, A ^ co. 



In other words, 7i:,„i„/ (v?) measures the minimal extent 
of the worldHne in the x direction of a coordinate sys- 
tem rotated by the angle (p. In Eq. ([25]), we also en- 
counter 7z„ax,f(^)) denoting the x-dependent envelope of 
the worldline in positive z direction. All these geometric 
properties of a worldline are displayed in Fig. [31 

The Qi.p.,ip functional in Eq. l[25|) generalizes the case 
of perpendicular plates (_L) for ip = ■k/2 and the case of 
one semi-infinite plate parallel to a infinite one (Isi) for 
if 0; both edge configurations were studied in detail 
in [ll,|4i|. 

Let us define 

7ni(x) EE XCOs{(p) + Sm{(p)'-fz^^^{x) - 7i^mi„(<<5)- (27) 

Inserting Oi.p.,,^ for f ^ into Eq. (flTI) leads to the 
Casimir energy density of the incHned plates 

E^-P-^'fi _ csc((^) 

^ ^ {in)D/^D-l){D -2) 

X /^J^^'"^^ dx ^^-\x)y (28) 

Equation ([28]) is shown as a function of ip in Fig. [H for 
Z? = 4. For (p = Tr/2 and D = 4, we rediscover the per- 
pendicular plates result (l6l . [45| as a special case. Inci- 
dentally, the integral in Eq. ([281 can be done analytically 
for = resulting in {X^/D) ^ {D~l)r{D/2)C{D). To- 
gether with the (^-dependent prefactor, Eq. ([28]) diverges 
as ^ as it should. This is because Eq. ([28l) corre- 
sponds to the energy per unit edge length, whereas for 
(p ^ the Casimir energy becomes proportional to the 
area of the semi-infinite plate. We devote the whole next 
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FIG. 4: Normalized Casimir energy per edge length — x 

a^~^ sin((p) of the inclined-plates (i.p.) configuration in D = 
4 versus the angle of inclination ip. For ip = Q, this function 
can be evaluated analytically, yielding 7r'^/2880 ~ 3.427- 10"'^. 
We have used 20000 worldlines with 10^ ppl each. 



section to analyzing how the Hmit ^ yielding the Isi 
configuration can be obtained. 



C. Inclined plates, </? — » limit 

It is instructive to study the Hmit of a semi-infinite 
plate parallel to an infinite plate (Isi), — ^ 0, as it 
involves a subtle limi ting process. Recalling the gen- 
eral considerations of [laSal for the Isi case, the total 
Casimir interaction energy decomposes into 



E. 



Isi, edge 



(29) 



where i?c^''"/A is the usual Casimir energy per unit area 
of two parallel plates Eq. (|23| . with A being now the 
area of the semi infinite plate. The so called edge energy 
^isi.odgc jiieasures the contribution that arises solely due 
to the presence of the edge. 

In the limit ^ 0, this decomposition is naturally 
achieved by inserting Oi.p.,<^=o of Eq. (|25l) into Eq. ifTT]) 
and performing the z integral first. This leads to 



^lsi,edgc 



1 



(47r)^/2(L>-2)a 



D-2 



(30) 




D-2 



This representation can straightforwardly be computed 
numerically Of course, for truly infinite plates, the 
edge effect being proportional to the length of the edge is 
completely negligible in comparison with i?c the lat- 
ter being proportional to the area of the plates. However, 
deahng with finite plates, the edge effect contributes to 
the Casimir force, effectively increasing the plate's area 

m. 



Of course, the same result has to arise from the general 
inclined-plates formula Eq. (|28|) in the Hmit (p ^ 0. How- 
ever, this representation naively exhibits a divergence in 
this limit. To find the origin of the divergence, we de- 
compose Eq. ((28|) into the parts corresponding to the 
edge effect E^'^^'^^'^ and the semi-infinite-plates energy 



i?i , characterized by the integrals /_ 



. dz and 



lo 



LJ2 



.dz, respectively . Here Lz/2 denotes the (infi- 



nite) length of the semi-infinite plate in z direction. The 
result for E^'^se,^ reads 



e: 



iedge,t/5 




'"'^"^^^ lz.^.A^h^'H^)dxy (31) 

which becomes Ef<i° in Eq. ((30l) as 0; Eq. ^ is 
therefore valid for = 0. 

On the other hand, if we naively expand the result for 
i?!'''' for small ip, we obtain 



e!1 



-{\^)Ly 



{4TT)D/2(D-2){D-l)DaD- 



'■Lp 



0(v5), (32) 



which is only valid for 7^ and does not reproduce 
Eq. (|22l) in the limit ^ 0. Instead of the energy per 
area, we have obtained the energy per length, which of 
course diverges in this limit. In order to rediscover the 
Casimir energy for the Isi configuration, the Hmits <p — > 
and the implicit limit Lz ^ 00 have to be taken in 
the right order. In Eq. ((32l) . the limit Lz ^ 00 has 
implicitly been performed first, which precisely leads to 
the divergence of the energy per edge length. Therefore, 
we need to first perform the Hmit ip ^ at finite Lz 
in order to obtain the desired energy per area. Starting 
from Ec'"^ at smaH (p, 



2(47r)^/2 



7x 



X / dx9 I -a 




dT 



7-(Z5+l)/2 



dz 



ZLp - 



we do the T integral first and obtain 



{iTT)D/^{D-2){D^l)Dp ' 



(33) 



(34) 



For small (Lzip), i.e., finite Lz and (/s 0, we can expand 
the last factor in ip, 



(477)^/2(1?- 1)L» 



1 1-i 



which for <p = corresponds exactly to the parallel- 
plates contribution i?c''''" in Eq. (|22l) . From Eq. l(34|) . we 
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FIG. 5: Normalized Casimir torque per edge length — x 

a^~^ of the inclined-plates configuration in D = 4 and its 
expansion around ip = n/2 and ip = 0, respectively, versus 
the angle of inclination ip. We have used 10000 worldlines 



with 10 ppl each. 



also observe that the other order of Umits, taking first 
Lz —>■ oo while keeping ip finite, reproduces the diver- 
gent behavior of the energy per edge length in Eq. l(32|) 
(as long as Re[D] > 2). The proper order of limits is 
similarly important at finite temperature with the addi- 
tional complication that another dimensionful parameter 
occurs. 



D. Casimir torque of inclined plates 

The Casimir torque DJ. P ''^ referring to rotations of one 
of the plates about the edge axis can easily be obtained 
by taking the derivative of the Casimir energy (|28)) . (or 
Eq. l(34| for small ip) , with respect to the angle of incli- 
nation: 



dip 



(36) 



For ip near tt/2, we can even set djx^i^i'p) / dip = be- 
fore taking the average with respect to the loop ensem- 
ble, simplifying the calculations. This is, because the 
derivative d'y x {ip) / dip changes its sign for perpendicu- 
lar plates = 7r/2 if the worldline is rotated by an angle 
TT about the normal axis of the lower plate, see Fig. [TJ 
Therefore, the sign correlates with the position of the 
minimum on the x axis of the lower plate. But the posi- 
tion X of the minimum 72-,„i„ {ip) does not correlate with 
the value of the integral in Eq. l(28|) leading to a mutual 
cancellation of terms involving d'-fx^i^{ip) / dip . 

For ip < 7r/2, we have to rotate the worldline about the 
normal axis of the inclined lower plate. Then, the cor- 
relation between the position of the minimum and the 
involved integrals does not vanish any more since the 
original and rotated worldline contribute differently to 



the integral. In general, expressions containing deriva- 
tives of 72;„i„('p) cannot be neglected even a,t ip = tt/2. 
Since the worldlines are not smooth, the convergence of 
averages of such expressions will be very slow. This is the 
case when calculating the coefficients of an expansion of 
Eq. (|36l) near ip = tt/2. Since the second derivative al- 
ready appears in the first expansion coefficient, more con- 
fident values are obtained by a numerical fit to Eq. (|36|) . 
There, only the first derivative is present. For Z? = 4, we 
obtain (see Fig. [5]) 



Li, 



0.00329 



0.0038 {^-ip^^ . 



(37) 



This should be compared to the worldline average based 
on the expansion of Eq. ((36l) around tt/2: the linear co- 
efficient in Eq. ll37| then yields 0.003 ± 0.0002. If we ne- 
glect all derivatives of 7a:„i„ {ip) the worldHne result reads 
0.00285 ± 0.00003. In all three cases 10000 worldHnes 
with 10^ ppl were used. 

For </? ^ 0, the Casimir torque diverges. The expan- 
sion about ip = {) can easily be obtained analytically from 
131 



£)i.p.,ip^O 



LyV{D/2)aD) 



{Att)I^/^{D -2)a 



D- 



2 2 ' 



(38) 



where we have used Eq. For = 4, Eq. ([Ml) yields 

Ly7rV2880aV^ ~ 0.00343Ly/aV^ being excellent ap- 
proximation to Eq. ([361 for ip not too close to tt/2. 

The divergent Casimir torque per length can be con- 
verted into finite torque per unit area by means of 
Eq. ([34|). Note that Eq. ([34|) leads to the classical re- 
sult for the torque. 



AL,T{D/2)(:{D){D^l) 
2{4TT)D/^aD 



(39) 



where A and Lz denote the semi-infinite plate's area and 
extent in z direction, respectively. For D = 4, Eq. ([39]) 
becomes AL^tt'^ /960a'^ « 0.0l03ALz/a^. 

A new characteristic contribution emerges from the 
edge effect Eq. ([31]). Unlike the total incHned-plate 
Casimir energy S^p-^ = si''^ + Ef'^'^''^, the edge en- 
ergy ([3T|) decreases with the angle of inclination ip, see 
Fig. [6l This leads to a contribution which works against 
the standard torque ([39]) . For D = 4, the correction to 
Eq. ([39]) emerging from the edge effect reads 



d: 



(0.003660 ±0.000038)^, 



(40) 



where we have used 10000 worldHnes with 10^ ppl each. 
The coefficient in Eq. (|40p was calculated by expanding 
Eq. ([3T|) around ip = 0. Equation ([40[) is shown in Fig. [6l 
We will see a similar subleading repulsive torque effect in 
the next section where we investigate finite-temperature 
contributions. 
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FIG. 6: Normalized edge energy per edge length — x 

of the inclined-plates configuration in D = 4 and its 
expansion around ip = versus the angle of inclination f. We 
have used 10000 worldlines with 10® ppl each. 



IV. FINITE TEMPERATURE 

Decomposing the Casimir energy at finite temperature 
T = into its zero-temperature part Ec{0) and finite- 
temperature correction AEc{T), 



E,{T)=E,{0) + AE,iT), 



(41) 



is straightforward in the worldUne picture by using the 
relation (flSl) . The finite-temperature correction is purely 
driven by the worldlines with nonzero winding num- 
ber. As the winding-number sum does not take direct 
infiuence on the worldline averaging, the complicated 
geometry-dependent part of the calculation remains the 
same for zero or finite temperature. This disentangles 
the technical complications arising from geometry on the 
one hand and temperature on the other hand in a conve- 
nient fashion. The same statement holds for the Casimir 
forceFc(r) =^^e(0) + AFe(r). 



A. Parallel plates 

In order to demonstrate the simplicity of the worldline 
method, let us calculate the well-known thermal contri- 
bution to the Casimir effect for parallel plates. In the 
following, we use the dimensionless parameter 



e = ar. 



(42) 



which distinguishes between the high-temperature ^ ^ 1 
and low-temperature ^ ^ 1 parameter region. 

Evaluating the general worldline formula for the 
Casimir energy Eq. lfT4|) using the parallel-plates Q func- 



tional of Eq. ^ results in {D > 2) 
£;ll(0) nD/2)C{D) 

+ \z^^riD/2)aD) 



2 2 



(43) 



4^2 



where the exponential integral function En{z) is given by 



En (z) = 



dt 



(44) 



and A again denotes the maximum extent of the world- 
line in the direction orthogon al to the plates. In the 
low-temperature limit, 2^ <$: (A^) = 7r/-\/3, the expo- 
nential integral functions vanish exponentially and can 
be neglected. We then obtain the small-temperature cor- 
rection to e\ (0) in D dimensions fully analytically: 



A^|(^-.o) _ r 



- D-l ■ 



^C(D - 1) (20 



eI{q) 



(45) 



The term (2^^ agrees with the standard textbook result 
It dominates the thermal correction to the Casimir 
force, yielding a comparatively suppressed power law de- 
pendence on the temperature, AFc(T) ~ for small T. 
This is an immediate consequence of the gap in the rele- 
vant part of the fiuctuation spectrum in this closed geom- 
etry. This term can also be understood as an excluded- 
volume effect: the volume in between the plates cannot 
be thermally populated by photons at low temperature 
due to the spectral gap. 

Incidentally, the leading contribution to the energy ~ 
is much less known. It does not contribute to the 
Casimir force, since it is independent of a when multiplied 
by the normalization prefactor Ec (0) . As we will see in 
section IIV B[ a-independent terms in the energy should 
not be viewed as mere calculational artefacts but can also 
contribute to observables such as the Casimir torque. To 
the best of our knowledge, Eq. l(45|) represents the first 
exact analytic formula for this leading small temperature 
correction to the free energy. 

The high-temperature limit of (|43l) can be obtained 
by a Poisson resummation of the winding-number sum 
(which is identical to returning to Matsubara frequency 
space). Our result agrees with [3] and reads: 



A£||(e^oo) 
e!{o) 



2r(i^)C(i^-i) V^ 



e (46) 



For arbitrary £, and D > 2, Eq. (|43l) can be evaluated 
numerically. Figure [7] shows the worldline result together 
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FIG. 7: Parallel plates: temperature dependence of the ther- 
mal contribution to the Casimir energy Ai?| (0) nor- 
malized to the zero-temperature result in D = 4-dimensional 
spacetime versus the dimensionless temperature variable ^ — 
aT. The worldline result for 1000 worldlines with 2 x lO'' 
points each is plotted together with the analytic expressions 
l|45|l - lf47ll . 



again revealing the power-law suppressed temperature 
dependence which is characteristic for a closed geome- 
try. Both magnitude and sign of the thermal force cor- 
rection can be understood as an excluded- volume effect: 
as the temperature is small compared to the spectral 
gap, thermal modes in-between the plates cannot be ex- 
cited. Hence, the thermal Stefan-Boltzmann energy den- 
sity outside the plates is not balanced by a thermal con- 
tribution inside. Thermal effects therefore enhance the 
attractive force between the plates. 

Let us finally remark that the comparison between the 
small-temperature limit of (|47|) (calculated with the help 
of the Poisson summation) and our analytic formula l(45|) 
closes a gap in the literature. With this comparison, 
we can find the exact value of the integral occurring in 
the prefactor of the leading low-temperature term in the 
energy. 



d.T^-^ [—2 -|- a;(cotha; + x csch x) 



C(3) 
2^2 ■■ 



(49) 



with the analytic asymptotics l(45|) - (|46|) and the known 
exact analytic formula for D = 4, see e.g. [i^ : 



coth(2n7r^) + 2n7r^ csch^(2n7rf) 



(47) 



Note that the result obtained in [46| for the electromag- 
netic field is twice as large as the result for the scalar 
field gll). 

Of course, taking the derivative of (|43|) . (|45l) or (06]) 
with respect to a also gives immediate access to the ther- 
mal corrections AFc(T) to the Casimir force. For in- 
stance, the low-temperature hmit results in 



Ai^c(aT< 1) 



Oa aT«l 

nD/2)aD)A 



(48) 



numerically corresponding to « 0.060897. This result has 
been observed numerically in the sum over odd refiection 
contributions to the parallel-plates Casimir energy in the 
optical approach to the Casimir effect 



B. Inclined Plates 



Whereas the inchned-plate geometry is much more dif- 
ficult to deal with than the parallel-plate case when using 
standard methods, there is comparatively little difference 
in the worldhne language. Inserting the inclined-plates 6 
functional l|25l) with 7ni [x) as in Eq. l(27|) into the general 
worldline formula (fT4l) yields 



Ai?i-P-'^(C), 



(50) 



where 



Ai?i-p-'^(o = - 



Ly csc((^) 
(47r)^/2a^-2 



(2e)^"^C(£'-2)r ' ^ ^ 



Ax^Ux))-aD-l){2if-^T ' ^ ^ 



E 



Ax 7^ '{x) X 

r 



2 2 



4^2 



4^2 



(51) 



is the thermal contribution to the energy. Here and in the D > 3 where all expressions exhibit well-controlled con- 
following, we confine ourselves to spacetime dimensions 
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vergence. In the low-temperature limit, the exponential 
integral functions can be neglected as long as 7111(2;) 7^ 
for all X. This is certainly the case for 7^ 0, but not 
necessarily for = 0. The latter case is again identical 
to the semi-infinite plate parallel to a infinite one, and is 
being considered separately in the next section and also 
in the appendix. For 7^ 0, the low-temperature limit is 
then given by the first two terms (first line) of Eq. (|5T|) . 
Note that the first ^gj-j^j does not contribute to the 

Casimir force, since it is an a-independent contribution 
to Ec if read together with the normahzation prefactor. 
From the second term, we obtain the low-temperature 
thermal correction to the Casimir force upon differentia- 
tion with respect to a. 



can be well approximated by 



i. p. ,1/37^0 



-Ly csc((p) 



(52) 



The temperature dependence differs from the parallel- 
plates case by one power of T, implying a significantly 
stronger temperature dependence at small temperatures. 
This is a direct consequence of the fact that we are deal- 
ing here with an open geometry. We emphasize that the 
result has been obtained fully analytically. In D = 4 and 
</? = 7r/2, our result agrees with the perpendicular-plates 
study of [l3| where this nontrivial interplay between tem- 
perature and geometry has been demonstrated for the 
first time. As shown therein, the thermal correction for 
this open geometry at experimentally-relevant large sep- 
arations can be an order of magnitude larger than for a 
closed geometry. 

Whereas the a-independent first term of Eq. ijSTj) does 
not contribute to the force, both terms in the first line of 
Eq. ((HJ contribute to the low-temperature limit of the 
Casimir torque. The thermal contribution to the torque 
is 



Ai5i-p-^(o = 



Alp 



(53) 



which at low temperature reads 

Ly cos(if) 



AL'i-P-'^(e^O) = - 



(47r)^/2 sin2((^) 



X [(7^„,n('/')Ax) -tan((^)(7^_^.J(/3)Ax)] 
+aaD - 1) {2Tf-' r V^) , (54) 



where J^"'"""' xdx^ ~ has been used, and we have 

introduced A^- = 7a;„^^ — 7a;„i„. This expression de- 
pends on only one nontrivial worldline average. In lim- 
iting cases, this average can be given analytically, as 
it reduces to the case described by Eq. (|24|) : we find 

-(7x™(¥' ^ 0)Ax) = (A2)/2 = 7rV6 and -(7.„.„('^ ^ 
'k/2)Xx) — (Xx)'^ /2 — 7r/2. For arbitrary cp, this average 



(7x,„,„(<^)Ax) 



--sin2(^) - — cos2((^), (55) 



as we will explain in the following. With Eq. (|26l) . we 
can write 

7rc„i„ (<p) = Ix (i) cos{ip) + {i) sin(<^) , (56) 

where t denotes the value of t that satisfies the minimum 
condition in Eq. (|26|) . Together with 



7;^.„.„(v) = -lx{i) sin{ip) + j.ii) cos(v?), (57) 

we can interpret {'yx,nin{f)i'lz,nin{f)) as the coordinates 
of the point {jx{t),Jz(t)) in the (^-rotated system. 

Since the jx and jz coordinates of each loop are gen- 
erated independently of each other, 7z(t) and Ax are not 
correlated. We therefore obtain 

(Ax7x„,i„(¥')> = (Aj;7x(t)) cos{ip) + {Xx){'Jz{i)) sm{ip). 

(58) 

By symmetry, the average (7z„i„('p)) vanishes, and we 
get from Eq. ([521) 



(7x (i)) sin((^) = (7, (£)) cos(^) . (59) 
On the other hand, (7x„,i„(v)) = —{^x)/'^- Substituting 



Eq. (|59|) into the average of Eq. (|56|) leads to {jzit)) = 
— sin((^)(Ax)/2 and {jxii)) = — cos{(p){Xx) /2, such that 
the desired Eq. (|55|) can be motivated by Eq. (|58|) . We 
would like to stress that only the second term in Eq. l(55|) 
has been estimated with the constraint imposed by the 
exactly known result for (p ^ 0, see above. The first term 
is exact and dictates the behavior of the perpendicular- 
plates limit. This result is compared to the numerically 
obtained data in Fig. [8l 

For the v-loop algorithm [l^l used here to generate the 
loops, the expectation value of the maximal extent A^^ is 
systematically smaller. This er ror is about the average 
spacing {\'rx{ti+i)-fx{U)\) w x/lf^N, see Eq. Asa 
consequence, the systematic error of (A^)/2 and (Ax)^/2 
is about (Xx)x/l/2N « 1.3 x lO^^ at TV = 10^. We 
observe a good agreement of the data with Eq. l|55p at 
(fi = and ip = tt/2 if the systematic error is taken into 
account. However, the agreement is actually perfect for 
all ip when using worldline estimates for the prefactors in 
Eq. l(55|) instead of 7r/2 and 7r^/6, see the modified curve 
in Fig. [HI This shows that Eq. (|55)) will well fit the data 
in the continuum limit N co. 

Let us return to the calculation of the Casimir torque. 
In the vicinity of the perpendicular-plates configuration, 
ip = 71/2— Sip, we can now obtain an expression to first 
order in Sip: 

LyScp 



^jji.p.,cp = TT/2-Scp _^ g-j 



(-C(^-2)r 



(47r)^/2 
D-2 



{2T) 



D-2 



lx„ 



2aTC{D - l)r 



D - 1 



(60) 
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FIG. 8: Angle dependence of the worldline average occur- 
ring in the Casimir torque l[54|l : {7a;„i„ (</?)Aa;). The numerical 
result is compared to the estimate l|55p (dashed line). The 
worldline result shows a small systematic error due to the fi- 
nite discretization. We can include the systematic error into 
l|55p by taking worldline estimates for the boundary values at 
if = and = n/2 (i.e. 1.567 sin^(v9) + 1.641 cos2(v5)), re- 
spectively, which are smaller than it/2 and 7r^/6. This yields 
the dot-dashed curve. The worldline result has been obtained 
from 5 10'' loops with 10^ ppl each. 



Here we have used (7^^ . (7r/2)A2^) — 0. Apart from 
{j'-^^.^{tt/2)Xx), Eq. l(60|) is an analytical expression. Us- 
ing'ISl), we obtain -(7;'^^j7r/2)A^) « ttVS-tt « 0.148, 
which is about ten percent of the dominating analyti- 
cal term in square brackets ^ 7r/2. We observe that the 
first term, which dominates in the limit aT 0, gives a 
contribution to the torque which drives the system away 
from the perpendicular-plates case ip = ■jt/2. Zero- and 
finite-temperature contributions thus have the same sign. 
The fact that = 7r/2 is a repulsive fixed point is also in 
agreement with naive expectations. 
For D = 4, Eq. ^ reads 



= dVT^ (0.0716 - 0.0957O , 



(61) 



which should be compared with the first-order term aris- 
ing from the T = contribution Eq. (jST)) . which reads 



i.p. ,93— TT /2—Sip 



obtain to first order in Sip 

^^i.p.,<p=V2-5v(^_^0) 



/Ly w 0.00329<5^/a2. Thus, for £» = 4 we 



FIG. 9: Integrals appearing in Eqs. I|5ip and II63II for 
D = 4. We have defined c^,To = ^X^^""*" da; 7ni(a;)^ and 

Cv,T^ = i^r^'^""""^ 7m(a;)^; see also Eqs. (HJ, (0) and 
Employing Eq. I|24p . we can evaluate c^=o,To = C(2) ~ 1.645 
and c^5=o,Too = \/^C(3) ~ 2.131 analytically. We have used 
10* worldlines with 10 ppl each. 



The high-temperature limit of Eq. ((511) can again be 
obtained by Poisson summation. The result is: 



A£:^-p-'^(^ ^ 00) = -£:i-p-'^(0) 



„D-2 



{x)(lx 



(4^)D/2aD-2(£,_3)(£)_2) sin((/j) 



(63) 



The remaining worldline average in this expression yields 
some positive finite number. Irrespective of its precise 
value for a given angle ip and D (the precise value of the 
integral for a specific (p can be read off, for instance, from 
Fig. [9] for either £» = 3 or £» = 4 and Fig. Hfor £> = 5), 
we stress that we observe the same linear dependence on 
temperature ^ = aT as in the parallel plate case. This is 
nothing but the familiar phenomenon of the dominance 
of the zeroth Matsubara mode at high temperatures, im- 
plying dimensional reduction, as discussed above. This 
mechanism is obviously geometry independent. Also the 
Casimir force remains attractive also for high tempera- 
tures. 



D. 



i.p..ip—7r/2-'Sip 



(0) 



(21 8 — 29 1^) (62) Semi-infinite plate parallel to an infinite plate 



In the validity regime of the low-temperature expansion, 
^ = aT <C 1, the positive first term is always dominant, 
hence perpendicular-plates case remains a repulsive fixed 
point. Most importantly, we would like to stress that the 
quadratic dependence of the torque on the temperature 
T^ '7^£'-2 ^j^g general case) for the inclined-plates 
configuration represents the strongest temperature de- 
pendence of all observables discussed in this article. 



A particularly interesting example for the geometry- 
temperature interplay is given by the semi-infinite plate 
parallel to the infinite plate (Isi configuration). In this 
case, the angle of inclination (p in Fig. [29] is zero. Analo- 
gously to Eq. l(29|) . the finite-temperature Casimir energy 
can be decomposed as 



Isi.cdgc 



(0 + £;^ii(e), 



(64) 
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where £;c"''ll(T) = Ec''''\o) + AEc^'^'Ht) correspond to 
the standard parallel-plate formulas as given in Eqs. l(23|) 
and l(43|) . with A being now the surface of the semi- 
infinite plate. Approaching the Isi limit of A£'c''''"(T) 
from the inclined-plates configuration in the limit (p 
is again a dehcate issue, as the proper order of Hm- 
its and Lz oo has to be accounted for, see 



Sect. IIII CI As the analysis is technically involved (but 
the outcome obvious), we defer it to the appendix. Let 
us here concentrate on the temperature-dependent edge 
contribution AE^^'-'"^se(^T) = E^^'-'"^se{T) - E^^''^°'^s<^{Q). 
We set = in Eq. JM)) and evaluate Eq. The 
result is (here and in the following, we confine ourselves 
to 13 > 3): 



(47r)^/2a^-2 



(20^"\(-D - 2)r 




Note that the first term, being the main contribution to 
the Casimir energy at small T, does not contribute to 
the Casimir force since it is a independent. Contrary to 
the case with ip 0, the exponential integral functions 
cannot be neglected in the low-temperature limit, since 
the argument of En{z) becomes zero at the lower bound 
of the integral for any ^ > 0. This results in a correction 
~ ^-D-i+a^ with a > 0, to the low-temperature limit of 
the first term. 

Here, however, we concentrate on the last term, as it 
gives rise to the thermal correction of the Casimir force. 
In the low-temperature limit, we find 



Isi.odgc 



(0 



(47r)^/2a^-i 



(66) 




D-2 



71=1 



For small ^, the main contribution to the integral comes 
from its lower bound as the exponential function rapidly 
decreases for large argument. At the lower bound, we 
can take the worldline average (7z„,ax (7x„i„ + x)) first, 
yielding a smooth function. We expand the latter in a 
power series. 



(72.„ax(7:r„„„ +x)) c„a 



(67) 



n=0 



where the exponents a„ do not necessarily have to be 
integers. Inserting (|67l) into (|66|) leads to 



-AE 



lsi,cdgc 



(C-o) 



(68) 



^^""-^r(^^±4^)C(i3 + a„-i), 



T1 = 
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FIG. 10 

{7zmax (72^min + s^)) obtaiued by 5 IC worldlines with 10' 
points each compared to the fit function of Eq. ifGOjl obtained 
on the interval x = [0,0.7]. The error bars have been plot- 
ted ten times larger. The observed small-a; power law 
with ai ~ 0.74 directly translates into a non-integer small- 
temperature behavior of the thermal edge contribution to the 

force, ARl^'-'^'^Se _ 



(7zmax(7a;min)) = Cq = 0, slucc for s, glvcu worldline 
{lx{t),'^z{t)) there exists a corresponding worldline in 
the ensemble with {'^x{t), — 7z(i))- We conclude that the 
coefficient of the term vanishes. We determine the 

higher coefficients c„ from computing (72^^^ (7a;„,i„ + 2;)) 
in the vicinity of x = by worldline numerics. Figure 
[10] depicts the form of (7z„,ax (72:mi„ + x)) near the lower 
bound x = 0. A global fit to this function including two 
coefficients ci , C2 is given by 



(7zn,ax(7a.„.„ +X)) 



0.9132 (x(1.500-x))"-^''^'^ 
1.234x0-7423 _ 0.6106x1-^423^ (gg) 



where we have neglected exponentially suppressed con- 
tributions. For the lowest-order term, we obtain 



where we have kept a2 — ai = 1 fixed. The resulting 
thermal correction to the force is shown in Fig. ifTTj) 
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FIG. 11: Thermal contribution to the Casimir edge force 
in the Isi configuration, — AFc°'''''^^°(^), plotted for a = 1 
and D = 4. WL: worldline result Eq. I|66p obtained us- 
ing 1000 loops with 10® ppl each. LTi, LT2: leading and 
next-to-leading low-temperature corrections 0.1098^^'^''^^ and 
0.1098C^-^^^^ - 0.131881C''■'^''^^ respectively, obtained from 
Eqs. (ISHl) and ([Ml, using 5 x 10'' loops 10*' ppl each. HT: 
high temperature limit obtained from Eq. I|70p . using 5 x IC* 
loops 10'' ppl each; a fit to the HT curve is provided by 
-5.24062(±0.0222)10-^ -h 1.591(±0.004138)10"^C- The inlay 
displays a magnified interval ^ = [0, 0.2]. 



for D = 4, where we compare the full numerical solu- 
tion with different orders of the expansion ((69l) and the 
high-temperature asymptotics, see below. As the low- 
temperature asymptotics is directly related to the lowest 
nonvanishing coefficient ai, we have also performed lo- 
cal fits to the function (7z„ax (7a:min + 2^)) the vicinity 
of a; = 0. Depending on the fit window, the leading ex- 
ponent can grow up to a\ ~ 0.8. (Of course, the fit 
window must be large enough to avoid that the worldline 
discretization becomes visible; otherwise, the exponent 
trivially but artificially approaches ai ^ 1 as the dis- 
cretized worldHne is a polygon on a microscopic scale). 

In any case, we conclude that the low-temperature 
regime of the Isi edge effect is well described by a non- 
integer power law, AF^i^'-^^^ge ^ j^D-i-fai ^ j^d-o.s^ 

where the fractional exponent arises from the geometry- 
temperature interplay in this open geometry. Of course, 
our numerical analysis cannot guarantee to yield the true 
asymptotic behavior in the limit ^ — > 0, but our data 
in the low-temperature domain 0.01 ^ C ^ 0-4 is well 
described by the non-integer scaling at next-to-leading 
order. 

Let us finally turn to the high-temperature Hmit of 
Eq. (|65)) which can again be obtained by Poisson sum- 
mation. The result for the edge energy reads 

^_gisi,cdgc(^ ^ 00) = -E]^^^'"^^^^)) (70) 
2^Ly ( /o^^ dx x^-3^,„,,Jx + 7,„,.j) 

(47r)^/2a^-2(i^-3) ^' 



where the worldline average is subject to numerical eval- 
uation. The resulting high-temperature Hmit of the 
Casimir force is shown in Fig. I|70p for D = 4. The 
high-temperature limit is again linear in T in accordance 
with general dimensional-reduction arguments. 



V. CONCLUSIONS 

In this work, we have provided further numerical as 
well as analytical evidence for the nontrivial interplay 
between geometry and temperature in the Casimir ef- 
fect. Whereas closed geometries such as the parallel- 
plates case exhibit a comparatively strong suppression 
of thermal corrections at low temperatures, open geome- 
tries such as the general inclined-plates geometry reveal a 
more pronounced temperature dependence in this regime. 
The terminology open and closed corresponds to the ab- 
sence or presence of a gap in the relevant part of the 
spectrum of fiuctuations which gives rise to the Casimir 
effect. In closed geometries, the spectral gap inhibits siz- 
able fiuctuations at temperatures below the scale set by 
the gap. By contrast, open geometries allow for sizable 
thermal fiuctuations at any value of the temperature. 

Concentrating on the incHned-plates geometry in D 
dimensions, the temperature dependence of the Casimir 
force can become stronger by one power in the tempera- 
ture parameter (implying thermal corrections which can 
be an order of magnitude larger than for a closed geom- 
etry). The inclined-plates geometry is particularly inter- 
esting as the limit of a semi-infinite plate parallel to an 
infinite plate (Isi configuration) is somewhat in-between 
open and closed geometries: the open part of the spec- 
trum only arises due to the edge of the semi-infinite plate. 
Interestingly, the resulting thermal correction numeri- 
cally shows a power-law temperature dependence with 
a non-integer exponent ^ j^D-a.a^ 

The strongest temperature dependence ~ ^0-2 -^^ ^j^^ 
low-temperature limit occurs for the Casimir torque of 
the inclined-plates configuration. This is, because it 
arises from the leading thermal correction of the inter- 
action energy which contributes to the torque but not to 
the Casimir force. 

Our results have been derived for the case of a fiuctu- 
ating scalar field obeying Dirichlet boundary conditions 
on the surfaces. Whereas this model system should not 
be considered as a quantitatively appropriate model for 
the real electromagnetic Casimir effect, our general con- 
clusions about the geometry-temperature interplay are 
not restricted to the Dirichlet scalar case. On the con- 
trary, all our arguments based on the presence or ab- 
sence of a spectral gap will also be vaHd for the elec- 
tromagnetic case. Whether or not the case of Neumann 
or electromagnetic boundary conditions leads to different 
power-law exponents for the temperature dependence of 
the "geothermal" phenomena remains an interesting ques- 
tion for future research. 

In view of the fact that most (strictly speaking all) 
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experiments are performed in open geometries, e.g., the 
sphere-plate geometry, at room temperature, an analysis 
of the geometry-temperature interplay of these experi- 
mentally relevant configurations is most pressing. 



same limit for the thermal correction to the energy. The 
decomposition of the Isi Casimir energy into bulk and 
edge contributions can also be performed for the thermal 
corrections. 



APPENDIX A: INCLINED PLATES, ^ - LIMIT Api-.^A^ - A Fedge,,,..^ , Aplh^/r^ 

AT FINITE TEMPERATURE - [I;) + 

We have analyzed the (jC — * behavior of incHned plates 
at zero temperature in Sect. IIII CI Here, we consider the where 



Ly csc((^) 



D-2 



(47r)^/2a^-2 



C{D - 2)(2e)^-2A,7..„.„(^) - C{D - l)r 



D - 1 



4e 



(a;cos(v3)-7:r.„,„(¥'))-7m {x)E3^ 



)} 



(Al) 



with 7m(a;) = xcos{ip) + sm{ip)jz,-m_^{x) — 7a;„,i„(<p) and 



(47r)^/2a^^2 



D-2 



C(^-2)(20 



D-2 



J2{ Tx^^~'ix)E,_o 



4^2 



(A2) 



Whereas AEf'^'^^'Pi^) remains finite, A£:i'^'^(0 shows a 
divergent behavior as (y9 — > 0. Let us therefore concen- 
trate on AEc''^{^), in order to isolate the source of the 
apparent divergence which is related to the order of lim- 
its of — > OD and (p ~* 0. In the case of inclined plates, 
{Lzip) is infinite for all </? 7^ 0, resulting in a 1/(^5 divergent 
energy density (energy per length) for (/? — * 0. Parallel 
plates, on the other hand have a finite energy density 



(energy per area) and {Lz^p) = 0. 

In the following, we show how to obtain an analytic 
transition from (Lzp) ^ 00 to {Lzp) for small </? 
by working with large but finite Lz, and taking Lz 00 
at the end of the calculation. The first limit results in 
a divergent energy density per unit edge length of the 
inclined plates as v? — > 0, 



En 



2a^-2^(47r)^/2 



{X')T[^]Q{D-2){20 



r(f )C(^)(20 



D 



4^2 

r 



D-2 



E^ 



4^2 



D - 1 
2 

2 2 y 4^2 



(A3) 
0(1). 



The second limit corresponds to the finite energy density 
of exact parallel plates l(43|) . 



In Eqs. I| All! Asp , the z integration was performed first. 
Let us now do the proper time integration first. The 6 
function l(33|). valid for small Lp, refiects itself in the lower 




(A4) 



15 



0.01 




FIG. 12: Qualitative behavior of the thermal Casimir con- 
tribution of the bulk in the inclined-plates case Eq. I|A7P 
(red line) and its small-(/> (dashed blue line) and large-(/f> 
(dot-dashed magenta line) limit, respectively. Note the di- 
vergent l/4> behavior of the large-(^ limit which corresponds 
to the divergent (as (p ^ Q) energy per edge length in the 
inclined-plates formulas. For this illustration, we have chosen 
D = 4, a — P = \x ~ ^, ignoring the worldline average in 
Eq. I|A7P for simplicity. 



The proper time integration yields for Re[D] > 1: 



exp 



dT 



2ia + z(p) 
xnj3 



D-l 



(A5) 



D - 1 



~ \4(a -I- zipy 



Inserting Eq. ([AS]) into Eq. i(A4l) leads to 



7r^A£;< 



l,v-»o 



(a,/3) 



E 

n=l 

aXx ( r 



r 



D K^'P' 



2 ' (2a+4 



D-l Kn'l3' 
2 ' 4a^ 



-r 



D-l Kn-P\ 
2 ' (2a+0)^ 



D-l 



D-l Kn^P' 



2 ' {2a+<ji)' 



2{n(3) 



D-l 



' 2 ' (2a+0)^ 



4(/)(n/3) 



D-2 



Hr(f)-r(f-ia^ 



D 



^r(V)C(^-i) 



2/3 



D-l 



ar(f)c(iJ) 

RD 



(A6) 



where (j) = L^^p and A = LyLz/2. One can show that 
the first three terms of Eq. I|A6P are of order 0{4>^). 
The forth term is clearly 0{<j)). The last line can be 
converted into the parallel-plates energy density l(43|) by 
neglecting 4> with respect to a and using the identity 
z^Ei^aiz) = r(a,z); the error is of order The 
first-order correction to the parallel-plates case is there- 



fore encoded in the fourth term. The second-order cor- 
rection is in the first three terms since the (j)'^ terms cancel 
each other in the remainder. In this limit ^ 1), all 
sums converge for Rc[D] > 2. 

Let us rearrange ljA6P so as to investigate the (j> —* oo 
case with (p being small but finite: 
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7r^Ag'i''^(a,/?) 
A 




D-2 
2 n2 



D-l Xjn^ffi 
2 ' 4a2 



40(n/3) 



D-2 



L>-1 



(2a+0)^ J 



2(n/3)^-i 

) 



2 ' (2a+^ 



4a2 



(A7) 



The large-(/) behavior of the first two terms can be ob- 
tained through Poisson summation^ and reads 

(4a + (/))(ci +C2(2fl + 0)) (C3 + C4(2a + (/))) 

(2a + (/.)^ (2a + (/))^-i ' ^ ' 

where ci, . . . , C4 are constants, the values of which are 
of no importance. We see that Eq. (|A8P vanishes for 
D > 2. For D > 3 the terms vanish even if multipHed 



by the infinite length L^. Remember that the inclined- 
plates formulae at finite temperature are valid for D > 3 
as well. 

In order to keep the remaining terms of Eq. I|A7P finite, 
we multiply both sides with the infinite length con- 
verting the vanishing Casimir energy per area into the 
finite energy per length. The Poisson summation of the 
second line of Eq. (|A7P results in 



D-l \>^I3 
2 ' (2a+0)2 



(2a+0)2 



A?r 



CP -2) 



4^(n/3)^-2 

aA,r(^)Cp- 



^r(#)CP) 



ip{2a + (j)) 



di{2a + (t>) 



D-2 



ip{2a + 4>) 



D-l 



dr^ + c?e(2a - 
tpi2a + (j)) 



D 



(A9) 



where di,...,dG are constants. These terms contain- 
ing di's vanish for D > 3 and (j) ^ oo. Applying the 
identity z°i?i_a(z) = r(a, z) to the last three terms in 
Eq. l|A7p . we rediscover the inclined-plates formula (jASP 
from Eqs. (|A7p . (|A9p valid for small angles 93. From 
Eqs. (|A8P and (|A9p . one can infer that the first correc- 
tion to Eq. ([All is of order 0(l/(/)^"3). 
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